Can a mass inversion save solar neutrino oscillations 
from the Los Alamos neutrino? 
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Abstract 

^ In the light of the —> v e neutrino oscillations which may have 

been observed at the LSND experiment we explore the consequences 
of two inverted mass schemes where solar neutrino oscillations occur 
between v e and v T . The favored LSND value Am 2 = 6eV 2 leads to 
m Ve ~ m VT ~ 2.5 eV and m v ~ so that cosmology can benefit from 
a recently proposed "cold plus hot dark matter" structure formation 
^2 scenario with two equal mass light neutrinos (Cz/ 2 DM). Solar neutrino 

oscillations {y e — ► v T ) can occur with one of the large mixing angle 
solutions so that a serious conflict with f3/3 decay Majorana mass lim- 
_C its is avoided without invoking Dirac masses. However, there is a 

problem with the SN 1987A signal because of resonant T7 e oscil- 
&h lations which are expected to cause far higher v e energies at the 1MB 

and Kamiokande II detectors than have been observed. A small value 
Am 2 = 0.5 eV 2 at LSND, which allows for a relatively large v e -v^ 
mixing angle without conflicting with the KARMEN and BNL-E776 
experiments, would indicate m Ve ~ m VT ~ 1.62 eV and m v ~ 1.77 eV. 
This scheme of Cz/ 3 DM maintains, and even may improve, the essen- 
tial cosmological model implications for large-scale structure, leav- 
ing no conflict with SN r-process nucleosynthesis. It may improve 
the discordance between the SN 1987A neutrino spectra inferred from 
Kamiokande II and 1MB. 
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The LSND collaboration (Los Alamos) has recently confirmed early rumors 
that they were seeing evidence for — > u e oscillations [1], probably corre- 
sponding to Am 2 6eV 2 and sin 2 29 0.005 [2, 3]. This interpretation may 
be marginally consistent with previous exclusion areas of the BNL experi- 
ments 734 of Ahrens et al. [4] and 776 of Borodovsky et al. [5], as well as with 
the KARMEN experiment [6]. The atmospheric neutrino anomaly may be 
explained by — > v T oscillations with Am 2 10~ 2 eV 2 and nearly maximal 
mixing [7]. Taken together, these hypotheses hint at a neutrino mass matrix 
with and m v „ 2.5 eV. Primack et al. have shown that 

the presence of two 2.5 eV neutrinos as a "hot component" of the cosmic 
dark matter ("CV 2 DM") may provide an elegant resolution of the structure 
formation crisis that bedevils a pure cold dark matter cosmology [3]. 

This explanation for results from two neutrino experiments seems at- 
tractive, but at least two significant, if not necessarily fatal, objections can 
be raised. Whether the atmospheric neutrino anomaly can indeed be ex- 
plained by oscillations remains controversial as the best-fit parameters of 
Kamiokande [7] seem inconsistent with the exclusion plots of the Frejus [8] 
and 1MB collaborations [9]. Moreover, there is no room for solar neutrino 
oscillations v e — > or v e — > v T which need Am 2 10~ 5 eV 2 (MSW solution) 
or Am 2 10~ 10 eV 2 (vacuum solution). Thus, one needs to invoke a new 
sterile neutrino in order to explain the solar neutrino problem by oscillations, 
an ad-hoc assumption that we find very unsatisfactory. 

Instead, one could give up on — > u T oscillations as an explanation for 
the atmospheric anomaly and invoke u e — > u T oscillations for the Sun. This 
implies 777/j/ ~ m VT and so, it would destroy the "cold plus hot dark matter" 
cosmology if 777/j/ ~ 777/j/ ~ and m v 2.5 eV. Also, it implies a mass 
inversion with m v > m VT . Of course, if one is willing to contemplate a mass 
inversion one may equally consider the reverse situation with 
2.5 eV and m v ~ 0. For cosmology, this scenario would be equivalent to that 
of Primack et al. [3] while the solar neutrino problem is solved by v e — > u T 
oscillations. 

An additional benefit of the mass inversion between u e and is that 
it avoids a conflict [10] with the recent interpretation that r-process nucle- 
osynthesis occurs in the neutrino-driven wind emanating from the surface 
of a nascent neutron star a few seconds after the core collapse in a type II 
supernova [11]. The neutron star emits neutrinos of all flavors with energies 
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for which various authors find values in the ranges [12] 

flO-12MeV for i/ e , 
(E v ) = I 14 - 17MeV for F e , (1) 
( 24 - 27 MeV for v^ T and v^ T . 

This hierarchy of energies is caused by the different opacities for the different 
species: they are emitted from different layers of the neutron star. The rela- 
tion {E Ve ) < (Ep e ) implies that (3 processes shift the wind to a neutron-rich 
phase necessary for the r-process. The LSND mixing parameters imply the 
occurrence of resonant oscillations near the neutron star which would swap 
v e ^ and would thus invert the energy hierarchy to {E Ve ) > (Ep e ), causing 
a proton-rich wind [13]. With inverted masses m Ve > m v resonant oscilla- 
tions occur among anti-neutrinos instead so that z7 e <->■ z7 M . The resulting 
energy hierarchy (E Ve ) < (E^) < (E-p e ) = (E VfL ) = (E Vt ) = (E Vt ) should be 
compatible with r-process nucleosynthesis. 

A 2.5 eV mass for v e could conflict with Majorana mass bounds from 
neutrinoless (3(3 decay. Currently, the best limit comes from the Heidelberg- 
Moscow Germanium experiment with an upper limit of m Ve ^ 0.7 eV [14]. A 
conflict is avoided if the v e mass is of Dirac type, a solution which requires 
one to postulate the existence of low-mass right-handed neutrinos. More 
interestingly, one may retain Majorana masses if one recalls that what is 
really constrained in (3(3 decay experiments is the quantity 
3 

( m "e) = '52*j\U e j\ 2 m j (2) 

where mj (j = 1, 2, 3) are the mass eigenstates, U e j their admixture to i> e , and 
Xj = ±1 a CP-phase. By assumption, m 2 ~ and mi ra 3 2.5 eV, so that 
a negative relative 1-3 CP phase allows for a cancellation if the 1-3 mixing 
angle is large. Thus, it is enough to require (m Ve ) = 2.5eV(cos 2 9 — sin 2 8) < 
0.7 eV where we used U e \ = cos 8 and U e 3 = sin 8 with 8 the 1-3 mixing 
angle which in our scenario is the one relevant for solar neutrino oscillations. 
Therefore, it is required that sin 2 28 > 0.92, i.e., almost maximum mixing. 
However, the uncertainty of the nuclear matrix elements entering the Majo- 
rana mass bound may be as large as a factor of 2 [14] whence the limit may 
be as weak as (m Ve ) ^ 1.4 eV so that only sin 2 28 ^ 0.69 would be needed. 

The solar neutrino problem is solved by matter-induced oscillations if 
Am 2 10" 5 eV 2 and either sin 2 28 0.007 or sin 2 28 0.6 (e.g. Hata and 
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Langacker [15]), assuming a standard solar model with physical input pa- 
rameters within their recognized experimental or systematic uncertainties. 
Therefore, our scenario favors the large-angle solution. It would cause a pro- 
nounced day-night effect in the Superkamiokande detector which hopefully 
will commence operation in the spring of 1996; within a few months of data 
taking this diurnal signal variation could be identified [16]. If this were the 
case one would be led to speculate that the discovery of neutrinoless (3(3 decay 
must be imminent. 

Another large-angle solution is provided by vacuum oscillations with 
10 10 eV 2 (e.g. Krastev and Petcov [17]). In this case the oscilla- 
tion length is of order the annual distance variation between Sun and Earth 
so that ultimately this solution can be identified by a semiannual variation 
of the measured fluxes. (Because the flux from the monochromatic beryl- 
lium line would show a particularly conspicuous time variation this effect 
could be observed well at BOREXINO which is the only forthcoming solar 
neutrino detector sensitive to beryllium neutrinos.) The vacuum solution is 
very attractive in the sense that it allows for maximum mixing so that the 
(3(3 decay limits will never cause a problem. In fact, one would view u e and 
u T as two Majorana components of one Dirac neutrino except for the small 
1-3 mass difference which presumably would have to be explained by some 
higher-order correction to whichever physical effect that causes the neutrinos 
to have masses in the first place. 

The atmospheric neutrino anomaly remains unexplained by oscillations 
in our scenario unless one invokes a sterile neutrino for this purpose. A far 
worse problem, however, is posed by the neutrino signal that was observed 
from SN 1987A in the 1MB and Kamiokande II water Cherenkov detectors. 
Our inverted mass hierarchy avoids a "level crossing" between u e and on 
the way out of the supernovae, and so v e <->■ exchange is avoided, a process 
that would have suppressed r-process nucleosynthesis. By the same token, 
however, a level crossing does occur between z7 e and z7 M , and so their spectra 
are expected to be swapped. Therefore, the detectors should have measured 
higher-energy z7 e 's than have been observed. 

A detailed maximum likelihood analysis of the SN 1987A neutrino signal 
has been performed by Loredo and Lamb [19]. Among ten different ways 
to parametrize the SN neutrino flux they favored a cooling model where 
neutrinos are emitted from a sphere at a fixed radius with a temperature 
which falls exponentially in time, beginning with an initial value T . For the 
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present discussion the temporal evolution of the neutrino flux is of minor 
interest; we are mostly concerned with the fluence (time-integrated flux) at 
Earth as a function of the F e energy. If one assumes that the instantaneous 
neutrino spectrum is given by a Fermi-Dirac distribution with a vanishing 
chemical potential, the exponential cooling model yields (Ep e ) = 2.36 T for 
the time-integrated spectrum. With the help of this conversion factor we 
show in Fig. 1 Loredo and Lamb's confidence contours for (E-p e ) and the 
inferred total amount of emitted F e energy. Typical SN calculations yield 
Etot(v e ) ~ 0.5xl0 53 erg, in good agreement with the best-fit value of Fig. 1. 
The inferred best-fit value of (Ep e ) = 10.5 MeV is somewhat low relative to 
the predictions cited in Eq. (1). 

In Fig. 2 we show the expected spectra at the two detectors. The solid 
lines correspond to the exponential cooling model with the best-fit values 
(Ep e ) = 10.5 MeV (T = 4.5 MeV) and E tot (u e ) = 0.5xl0 53 erg. It agrees 
well with the dotted line which was derived from the numerical cooling cal- 
culation of Burrows [20], his model 55, for which the equation of state and 
matter accretion rate were adjusted such as to optimize agreement with the 
detector signals. In this model (Ep e ) = 11.0 MeV; the instantaneous neutrino 
spectrum was taken to be black-body. This comparison confirms that Loredo 
and Lamb's parametrization is a reasonable representation of what one might 
expect from a numerical cooling calculation, at least for the time-integrated 
spectrum. 

So far it was assumed that the instantaneous neutrino spectrum is a 
Fermi-Dirac distribution with vanishing chemical potential. However, the 
spectrum is probably "pinched" (suppressed low- and high-energy tails) as 
shown, e.g., by the Monte Carlo transport calculations of Janka and Hille- 
brandt [21]. These authors [22] have performed a maximum likelihood anal- 
ysis of the SN 1987A signal on the basis of a spectral parametrization of 
the fluence in terms of an effective temperature and an effective degeneracy 
parameter r\ where r\ > corresponds to pinching (suppressed high energy 
tail), 7] < to anti-pinching. For each detector separately, they found a 
best-fit value of r\ = 0; probably they would have found a negative value 
had they not imposed the restriction 77 > 0. For r\ = their best-fit re- 
sults are: (E-p e ) = 9.1 MeV and E tot (u e ) = 0.64xl0 53 erg (Kamiokande II) 
and (E^ e ) = 14.2 MeV and E tot (u e ) = 0.34xl0 53 erg (1MB). The confidence 
contours are not very restrictive on r\. Fixing r\ = 3.5 on the basis of 
what is theoretically preferred, the best-fit results are: (E^ e ) = 10.9 MeV 
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and E tot (u e ) = 0.53 xlO 53 erg (Kamiokande II) and (E-p e ) = 18.6 MeV and 
^tot(^e) = 0.23 XlO 53 erg (1MB). Even for rj = the best-fit (E-p e ) disagree 
substantially between Kamiokande II and 1MB; their 95% confidence con- 
tours barely overlap. This discrepancy is worse for the assumed pinching 
with 7] = 3.5. The data prefer anti-pinching in each detector separately, 
and to explain the relatively large number of events at 1MB (which is only 
sensitive to the high-energy tail) relative to Kamiokande II. Of course, anti- 
pinching of the time-integrated spectrum is not necessarily in conflict with 
the assumption of pinched instantaneous fluxes. 

The spectral distribution at the detectors lead to values for (Ep e ) which 
are slightly low, but compatible with predictions of numerical calculations. 
However, at a reasonable confidence level they are not compatible with 
the assumption that primarily converted F^'s had been observed. For F^'s 
the expected spectral pinching is far less pronounced than for F e 's so that 
the Loredo and Lamb analysis is essentially adequate. Their joint analy- 
sis between both detectors allows (E^ e ) = 15.5 MeV (T = 6.5 MeV) with 
^tot(^e) = 0.2 XlO 53 erg just barely at the 95% CL. From Fig. 2 (dashed line) 
the disagreement with the Kamiokande II signal is plainly visible. Even this 
large value for (Ep e ) is far below the expectation for converted F^'s. There- 
fore, the SN 1987A signal favors the normal over the inverted mass hierarchy 
by a large margin. 

The predictions of the spectral features of the neutrino signal are fraught 
with uncertainties. One is related to the neutrino opacities at high den- 
sity where one expects that fast spin-fluctuations lead to a suppression of 
the axial- vector scattering rate at densities beyond, say, 10% nuclear [23]. 
The magnitude of this effect has not yet been calculated from first princi- 
ples. Instead, the neutrino opacities have recently been "calibrated" from 
the SN 1987A signal [24]. Reduced neutrino opacities lead to a stiffening of 
the observable F e spectrum [24] and so, the average energies of Eq. (1) which 
were predicted on the basis of naive unsuppressed scattering cross sections 
likely are lower limits. Therefore, it is even less believable that converted 
F^'s were observed from SN 1987A. 

Does one really expect a complete swap F e <->■ F M by matter-induced os- 
cillations? For our inverted mass scheme the F e becomes degenerate with 
the F M when Am 2 /2E U = \[2G^n e with the Fermi constant Gp and the elec- 
tron density n e . With the LSND value Am 2 = 6eV 2 , with E v = 30 MeV 
which is characteristic of the detected neutrinos, and with an electron frac- 
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tion of 0.5 per baryon one finds that the cross-over occurs at a density of 
2.6 X 10 6 g cm -3 . This is far outside the neutrino sphere which is at a den- 
sity of around 10 13 gcm~ 3 . The oscillations are adiabatic if on resonance 
|Vlnn e | > (Am 2 /2E 1/ ) sin 26 tan 26. With the LSND value sin 2 26 = 0.005 
and E v = 30MeV, this condition reads |Vlnn e | _1 ^> 0.4km, which is easily 
satisfied even after the shock wave has expelled the stellar mantle so that the 
density falls off relatively rapidly near the neutron star surface. 

Of course, at relatively late times the resonance occurs at a few 10 km 
from the neutron star surface while at early times, before the shock wave 
has escaped, it occurs at a few 1000 km. Clearly, there will be complicated 
dynamics and perhaps large-scale motions of material and turbulence which 
render a simple spherically symmetric density profile to be too simplistic. 
While it appears implausible that the adiabatic condition will be destroyed 
for a substantial region of space or a significant period of time during the 
Kelvin-Helmholtz cooling phase of the neutron star, a more detailed analysis 
may well be called for before passing final judgement. 

Besides the electron density profile, one has to worry about neutrino 
self-interactions with regard to the matter-induced oscillation problem. The 
same issue arises for the normal mass hierarchy and the resulting v e <->■ 
oscillations, a problem that was studied in Ref. [13] in the context of a possible 
suppression of the r-process nucleosynthesis. The neutrino self-interactions 
seem to cause only a mild modification of the adiabatic conversion process 
[13]. Still, these issues have only begun to be fully understood, and so one 
may not be entirely satisfied yet with the approximate treatments that are 
hitherto available. 

The interpretation that the SN 1987A signal was indeed caused by z7 e 's 
by virtue of the u e p — > ne + reaction is in itself subject to doubt because in- 
stead of the expected isotropic signature one found a forward peaked signal, 
especially for the high-energy events, which is compatible with isotropy at 
far less than a 1% probability [25]. Because there is no plausible alternative 
explanation available at the present time one must settle for the interpreta- 
tion of a rare statistical fluctuation of the angular distribution of the events, 
a fluctuation which should not affect the interpretation of the average en- 
ergies and number of events. Still, one remains uneasy about far-reaching 
conclusions which rely on a signal with rather extreme anomalous features. 

In spite of these systematic uncertainties associated with the interpreta- 
tion of the SN 1987A neutrino signal it is difficult to simply discard it in 
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favor of our inverted mass scenario. Is there another way out? One possi- 
bility is indicated by the observation that the LSND experiment alone does 
not allow one to fix Am 2 with certainty. The £/E distribution of the events 
which are associated with the F M — > u e channel also allow for the values 0.5, 
2, 6, or 10 eV 2 with a reasonable significance [2]. The value 6eV 2 was favored 
because of a sensitivity dip of KARMEN and BNL-E776 [2, 3]. However, if 
the LSND-favored sin 2 28 is indeed below 0.01, contrary to early discussions, 
there may be no conflict with KARMEN at any Am 2 [6] . This leaves a poten- 
tial compatibility problem with BNL-E776 for which two different exclusion 
plots are in circulation. One 1 shows a sensitivity minimum at Am 2 = 6eV 2 ; 
it requires sin 2 28 < 0.005 at 90% CL. In contrast, the exclusion plot in the of- 
ficial BNL-E776 publication [5] shows a sensitivity minimum at 8 eV 2 (a value 
disfavored by LSND) where sin 2 28 < 0.003 is required at 90% CL. It shows 
high sensitivity at 2, 6, and 10 eV 2 where sin 2 28 ^$ 0.002. This nominally 
excludes most of the LSND-favored values unless the errors of BNL-E776 
were larger than stated, or if one appeals to a rare statistical fluctuation. 
On the other hand, if the LSND results were explained by Am 2 0.5 eV 2 
there does not seem to be even a nominal disagreement with KARMEN or 
BNL-E776 even for mixing angles as large as sin 2 28 0.02. It remains to be 
seen if the putative evidence for oscillations in the — > u e channel at LSND 
selects any of the Am 2 values with superior significance. 

On the face of it, a low Am 2 destroys the beautiful cosmological scenario 
of Primack et al. [3]. However, it is only required that the sum of the neu- 
trino masses is about 5eV. Together with Ira 2 — ra 2 I 0.5 eV 2 one could 
contemplate a mass matrix 1.62 eV and m v „ 1.77 eV. This 

CV 3 DM model has an interesting implication for cosmology: it provides a 
dip in the ratio of the power spectrum relative to that for the 1 neutrino 
model, at fixed Q, v , at approximately twice the scale found in a similar com- 
parison for the 2-neutrino model, that is, near 20/?. -1 Mpc, by analogy with 
the computations cited in reference [3]. This would help further suppress the 
abundance of massive galaxy clusters and amplitudes of large-scale flows, and 
thereby lead to better consistency with the observational constraints [27]. 

1 This plot has appeared in the Proceedings of Neutrino 90 (pg. 148), in the review [26], 
in the KARMEN paper [6], and was recently shown at the Moriond conference in Villars, 
Switzerland, January 1995. We understand that this plot was based on a narrow-band 
beam, while the 1992 publication on a wide-band beam (W. Lee, communicated to us by 
D. Caldwell). Therefore, the 1992 version [5] should be used. 
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Even with the proposed normal z^-z^-hierarchy the r-process nucleosyn- 
thesis is left unscathed as the resonant oscillations would now occur too far 
away from the neutron star surface; even the value Am 2 6 eV 2 was close 
to the edge of the sensitive regime [13]. 

An explanation of the solar neutrino problem by v e — > u T oscillations still 
requires a large mixing angle in order to avoid a conflict with the (3(3 Majorana 
mass limits. However, with 1.62 eV the nominal Germanium 

bound (m Ve ) < 0.7eV translates into sin 2 28 > 0.8 so that both large-angle 
solutions are possible without the need to appeal to large uncertainties of the 
nuclear matrix elements. 

With m Ve < m v there is no problem with the SN 1987A signal. On the 
contrary, if the solar neutrino problem is solved by the large-angle MSW solu- 
tion the u e -u T mixing parameters are such that matter-induced oscillations in 
the Earth are important; this effect leads to the diurnal signal variation that 
would be expected at Superkamiokande. For SN 1987A, the detectors would 
have observed a certain combination of primary z7 e 's with z7 T 's whose spectra 
would have been partially swapped by the Earth effect. As Kamiokande II 
looked at SN 1987A through approx. 3900 km of the Earth, 1MB through 
approx. 8400 km, the effective spectrum observed by 1MB would have been 
stiffened relative to Kamiokande II, partially explaining the "tension" be- 
tween the best-fit (Ep e ) inferred from the two detectors [28]. In our first 
mass scheme with m Ve > m v this effect would not change the relative spec- 
tra because the z7 e 's arriving at Earth would have been converted F^'s with 
a spectrum identical to that of the z7 T 's. 

In summary, if the LSND neutrino survives further experimental scrutiny, 
an explanation of the solar neutrino problem by oscillations among sequen- 
tial neutrinos requires a near u e -u T degeneracy, and so requires either a mass 
inversion between v e and or between and v T . If the LSND signature is 
explained by Am 2 = 6 eV 2 as originally proposed, the r-process nucleosyn- 
thesis argument as well as the CV 2 DM cosmology favor 777/j/ ~ 777/j/ ~ Jim 5eV 
and m v ~ 0, a scenario that is problematic in view of the SN 1987A neu- 
trino observations. If the LSND value is Am 2 = 0.5 eV 2 , cosmology favors 
1.62 eV and m v 1.77 eV which avoids any conflict with 
SN physics. Unless one appeals to neutrino Dirac masses, the (3(3 Majorana 
mass limits require in either scenario a large-angle solution to the solar neu- 
trino problem which likely can be identified at the forthcoming solar neutrino 
experiments by a diurnal or semiannual time variation. 
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Figure 1: Projection of the confidence volume in a maximum likelihood 
analysis of the SN 1987A neutrino signal, assuming an exponential cool- 
ing model (adapted from Loredo and Lamb [19] with our conversion factor 
{Ev e ) = 2-36 T ). 
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Figure 2: Expected and measured distributions of the events in the 
Kamiokande II and 1MB detectors from SN 1987A. On the vertical axis, 
1 event in the binned spectra corresponds to 0.2. The solid and dashed lines 
are predictions for exponential cooling models. For (Ep e ) = 10.5 MeV, which 
corresponds to the best-fit value of Fig. 1, a total energy of 0.5xl0 53 erg 
was assumed to have been emitted in z7 e 's. For (E^ e ) = 15.5 MeV we chose 
-E'tot(i'e) = 0.2xl0 53 erg, corresponding to the right-hand tip of the 95% con- 
tour in Fig. 1. The dotted line is the prediction from the model 55 calculation 
of Burrows [20]. 
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